Multi-speckle diffusing wave spectroscopy with a single mode detection scheme 
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We present a detection scheme for diffusing wave spectroscopy (DWS) based on a two cell geometry that 
allows efficient ensemble averaging. This is achieved by putting a fast rotating diffuser in the optical path 
between laser and sample. We show that the recorded (multi-speckle) correlation echoes provide an ensemble 
averaged signal that does not require additional time averaging. Furthermore, combined with traditional two- 
cell DWS, the full intensity autocorrelation function can be measured with a single experimental setup. The new 
scheme provides access to a large range of correlation times thus opening a new experimental window for the 
study of slowly relaxing and arrested systems, such as viscoelastic complex fluids, colloidal glasses and gels. 



The surging interest in slowly relaxing and arrested col- 
loidal systems such as gels or glasses 01 01 has created 
a need to monitor dynamic properties on time scales of sec- 
onds and minutes. Light scattering is certainly one of the 
best methods for this purpose since it offers convenient access 
to such key dynamic properties as the intermediate scattering 
function or the particle mean square displacement. Tradition- 
ally a single speckle mode of scattered light is detected and 
fluctuations are recorded over a time much longer that the re- 
laxation time. However, this time averaging scheme is not 
applicable to rigid, nonergodic systems. For these systems 
the ensemble average can be obtained by summing a collec- 
tion of consecutive experiments conducted on different sam- 
ple realizations. Usually the sample is translated or rotated 
and a scan over a large number of independent speckles is 
performed A major drawback of this approach 

is the extensive duration of measurements. It is not unusual 
today to investigate relaxation process on time scales of sec- 
onds and minutes with a corresponding measurement time of 
hours and days. As a matter of fact several authors have re- 
ported data collection time more than a day for a single in- 
tensity correlation function (e.g. |8|, |9[). Besides being te- 
dious and time consuming this approach is restricted to the 
systems in (quasi-)equilibrium. Only the advent of multi- 
speckle detection schemes made it possible to conveniently 
monitor very slow relaxation processes. Dynamic light scat- 
tering using a digital (CCD/CMOS) camera as a detector of- 
fers the possibility to perform simultaneously a large number 
of independent experiments thus achieving ensemble averages 
in real time 1 10]. Compared to single photon counting the lim- 
ited dynamic range of digital cameras and the typically high 
dark counts result in a somewhat lower accuracy. But unfor- 
tunately with a time resolution of typically 1-10 ms digital 
camera based detection is restricted to rather long correlation 
times. Thus, traditional photon correlation spectroscopy has 
to be made as well in a separate experiment if access to the 
full range of correlation times is required. 

All considerations above apply equally to both dynamic 
light scattering (DLS) in the single scattering regime and dif- 
fusing wave spectroscopy (DWS) in the multiple scattering 
regime llllfl^fl3ll . However, due to the strong multiple scat- 
tering DWS offers more flexibility in the experimental design, 
which we have exploited in our approach. In this paper we 
report on a new two-cell detection scheme for diffusing wave 
spectroscopy (DWS) that provides an effective multi-speckle 



averaging using single mode detection. To obtain an ensemble 
averaged signal we illuminate our sample with the laser light 
scattered from a rotating diffuser and we analyze reflected or 
transmitted light. We show that echoes in the recorded cor- 
relation function appear at any revolution while the correla- 
tion function of the sample remains finite. Each echo signal is 
generated by a large number of independent speckles thus effi- 
cient ensemble averaging is performed. Moreover, we demon- 
strate that the intensity correlation function of the sample can 
be extracted from the two-cell echoes. 

The detection of single and multiple scatterin g co rrelation 
echoes was discussed in previous articles Jl4lll5Lll'rlfl7ifl8|] . 
Echo DWS was initially introduced in the analysis of non- 
linear shear deformation Il4 1 1 511 . In the single scattering 
regime Pham et al. recently demonstrated the use of echo 
DLS for efficient ensemble averaging 0. Our DWS echo- 
scheme implements a new physical optical principle to record 
the ensemble averaged intensity correlation function. In con- 
trast to previous echo experiments in our case the sample is at 
rest. Thus high rotation or oscillation frequencies can be real- 
ized without any mechanical disturbance of the system under 
study. The possibility to perform multi-speckle experiments 
with a traditional DWS light scattering scheme thus opens the 
pathway for a new type of fast and precise experiments. If 
combined with the well established two-cell DWS technique 
(TCDWS) 1 20] correlation times from 10 ns or less up to dura- 
tion of measurement can be accessed. Such improved experi- 
mental performance is mandatory if progress shall be made in 
the expanding field of slowly relaxing and arrested systems, 
such as viscoelastic complex fluids, colloidal glasses and gels 
flii. 

Our experimental setup is shown in Fig. ^ A frequency- 
doubled Nd:YVo4 laser ("Verdi" from Coherent) operating 
at Aq = 532 nm is used to illuminate a circular ground 
glass mounted on a 5-phases stepper motor (RK-564 AC from 
VEXTA). Through scattering and dephasing the ground glass 
creates a speckle with a nearly gaussian optical field 1 12D . We 
collect the transmitted light coming from the ground glass and 
focus it onto the sample with a spot size diameter of roughly 
5 mm. The scattered light is then collected with a mono mode 
fiber and analyzed by a photomultiplier and a digital pho- 
ton counter (correlator.com, New Jersey, USA). The photon 
counting device records the time intervals between photons 
arrivals in a data file with a resolution of 1/60 MHz = 12.5 ns. 
The high temporal resolution, comparable to the dead time of 
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FIG. 1: Experimental setup: Laser light is scattered from a ground 
glass rotated by a fast stepper motor (1) and the transmitted light is 
collected by a lens (2) to illuminate the sample (3). Single mode 
fibers collect the scattered light either in transmission (4a) or reflec- 
tion (4b). The collected light is subsequently analyzed by a single 
photon detector and digital photon counter (5) 



the detector, assures that there is no more than one photon ar- 
riving at a time step for a typical experimental count rate of 
10-500 kHz. The stepper motor is operated at frequencies up 
to 75 Hz. The measured intensity correlation function (ICF) 
[g M (r) - 1] contains information from both the sample and 
the rotating diffuser. To distinguish the dynamics due to the 
sample internal motion [gf^ij) - 1] and due to the rotating 

(2) 

diffuser [g s (t) - 1] in a quantitative way we take advantage 
of previous studies of similar two-cell geometries. Scheffold 

(2) 

et al. have shown that the ICF g M (r) - 1 from a sandwich of 
two optically independent cells can simply be expressed by a 
product of the correlation functions of the two individual cells 
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The two cell geometry can be also realized using a very 
slowly rotating diffuser as suggested by Viasnoff et al. l2lll . 
In this version the two cells are separated by a distance of sev- 
eral centimeters which ensures complete decoupling of light 
propagation in both cells. This realization of Two-Cell DWS 
is similar to our experimental setup. However, in previous 
studies the diffuser was rotated slowly in order to average a 
large number of arrested speckles over time whereas in our 
case the diffuser motion is fast and periodical. To illustrate 
the different contributions we first consider the scattering sig- 
nal from a rigid sample without internal motion. Any fluc- 
tuation of the detected laser light is produced in this case by 
the motion of the random diffuser. Fig. Eta) shows the typi- 
cal intensity correlation function of light scattered by a teflon 
slab of 2 mm thickness whith a motor rotation frequency f r 
close to 40 Hz. At short times the rotation gives rise to a com- 
plete decay of the ICF on a characteristic time x r set by f r 
and a corresponding echo width of 2r r . Each of the speckles 
reappear identical in the next revolutions resulting in echoes 
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FIG. 2: (a) and 1st order correlation echo in backscattering from 
a rigid teflon slab for a rotation frequency of f r » 40 Hz. The echo 
width at 1/e points is 2r r « 1.15 [AS. (b) DWS echoes in backscatter- 
ing for a sample of TiOi in glycerol. Inset shows the echo shape for 
this sample. 



(2) 

in [g E (t) — 1]. At T r — l/f r and any multiple integer val- 
ues n = 2, 3,4, ... a correlation peak is observed. Echoes are 
found indistinguishable for backscattering and transmission 
(data not shown). A detailed discussion of the echo shape is 
beyond the scope of this letter but we expect many similar- 
ities to the formalism developed for single scattering echoes 
However, the number of speckles sweeping over the de- 
tector can be readily estimated from our experiments to be 
N ~ T r ]2r r > 2 ■ 10 4 . This means already after n + 1 revo- 
lutions the correlation function of the nth echo is known with 
an accuracy better than 1 / y/N ~ 1%. 

In the following sections we will discuss the application 
of the two cell echo and compare the different data analysis 
schemes. We have prepared a colloidal dispersion of titanium 
dioxide powder (Ref: 0041255 from Warner Jenkinson Eu- 
rope Ltd.), particle diameter roughly 200-300 nm, in glycerol 
at a volume fraction of m 0.5%. To further increase the viscos- 
ity the sample is kept at 5.7 ± 0.5° C . In this viscous opaque 
medium the correlation function decays over a range of lag 
times accessible both to the echo technique and traditional 
time averaging. Fig. [2] shows the result of the echo measure- 
ment in backscattering geometry for a frequency of 40 Hz. 

In the viscous glycerol solution the scatters undergo diffu- 
sive motion expressed by the mean square particle displace- 
ment (r 2 (r)) = 6Dqt, where Do = kgl ' j(mr\a is the Einstein 
diffusion coefficient. For the backscattering geometry the ex- 
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perimental ICF is described by the expression 1 12]: 

gf(T) - 1 = fi ■ exp [-2y V6t/t J (2) 

with a relaxation time to = 1 /DkZ and a factor fi that describes 
the intercept of the correlation function, y is a constant of 
order y « 2.1 - 2.3 (VH detection) lEl . 

It follows from Eq. Q that gf\r) - 1 can be extracted from 
the echo-peak value since at lag times r = «//,• one expects 
g E (n/fr) -1 = 1. However, the echo peak hight might be 
affected by slight imperfections in the rotation. Furthermore 
a detailed resolution of the peak maximum can be costly in 
computation time in particular for higher rotation frequencies. 
A more practical way of dealing with this problem is to ana- 
lyze the peak area rather than the peak height. Pham et.al 
have shown for the case of single scattering echoes that the 
peak area is directly proportional to ideal peak height, suffer- 
ing very little from slight imperfections in rotation I17I1 . The 
peak aera can be obtained by numerical integration or sim- 
ply by increasing the sampling time t s . The latter approach is 
equivalent to triangular-weighted integration |23] and more- 
over significantly reduces the computation time of g M (r) - 1 . 

The relevant information is contained in the sample corre- 
lation function g s (f) — 1. To study the influence of the inte- 
gration time window and the sampling time r s we have varied 
both parameters over large range. Instead of numerical inte- 
gration we simply sum the points of the correlation function 
g M (r) - 1. No significant dependence is found (within rea- 
son) as long as the time window covers well the correlation 
peak. As a matter of fact a single correlation channel of sam- 
pling time t s = 12 \is centered at the echo position provides 
nearly the same level of accuracy as integration over 960 chan- 
nel at 12.5 ns resolution. In Fig.[3]we compare the echo-data 
collected during 12 seconds to a time averaged measurement 
over 20 minutes. Dividing by the normalization factor fi we 
obtain basically identical results from both methods|26J. De- 
spite a dramatically shorter measurement time, the noise level 
is lower for the Echo-measurements. 

In order to record these narrow correlation echoes special 
emphasis has been given to the data analysis. In a traditional 
linear correlator the ICF is obtained by averaging products 
of photon counts in a certain sampling time interval t s sepa- 
rated by lag time r (channel) 1 1JJ. Over several decades of 
lag times the linear channel layout however becomes imprac- 
tical and the computation time increases rapidly. The usual 
solution to this problem in hardware correlators is to use a 
multi-tau scheme as introduced by Schatzel 1 23]: to decrease 
amount of data the sampling time is doubled after a fixed num- 
ber of steps. The correlation function is well resolved for the 
shortest lag times while the resolution is decreased for longer 
lag times. This loss in resolution is however not acceptable 
for our signal since all echoes have the same width, even at 
large lag times. On the other hand only a small number of 
correlation channels is needed to resolve the ICF in the vicin- 
ity of the correlation echoes. Furthermore not all echoes need 
to be resolved to cover a given range of lag times. Thus, if the 
parameters are optimized, as shown later in the text, simple 



multiplication can still be the most efficient method. 

There are alternative techniques that allow a fast compu- 
tation of the correlation function. One of the most efficient 
methods is based on the fast Fourier transform (FFT). Ac- 
cording to the Wiener-Khinchin theorem the auto-correlation 
function of the signal can be determined as an inverse Fourier 
transform of its power spectrum which is the square of the 
Fourier amplitudes (for details see 1 24]). For a given sampling 
time t 5 FFT provides the correlation function for all available 
channels much faster than any other technique. However, in 
order to to resolve the echoes shape the sampling time has to 
be decreased which requires large data arrays and long pro- 
cessing time. 

Quite a different approach to calculate the ICF was put for- 
ward by Chopra and Mandel |25] in their photon time-of- 
arrival correlator. It estimates the distribution of time intervals 
between photons arrivals which is shown to be proportional to 
the ICF. In contrast to the methods mentioned above compu- 
tation can not be accelerated by increasing the sampling time 
since photons arrivals have to be recorded as separate events. 
This approach is much faster only if high resolution data is 
required and if the data is provided in a suitable format as it is 
the case for our photon recorder. We have used this method to 
calculate the ICFs in Fig.[2]at the highest available resolution 
of 12.5 ns. 

To optimize the data analysis we compare different meth- 
ods: the linear correlator, FFT and the time-of-arrival correla- 
tor. The data of a 12 seconds measurement at an average count 
rate of 210 kHz is analyzed under identical conditions. Cor- 
relation echoes are calculated up to echo number 240, corre- 
sponding to r = 6 seconds and the time delay between echoes 
is doubled after each linear block of 16 and a total number of 
60 echoes is calculated. The time resolution of time-of-arrival 
correlator is determined by the 12.5 ns hardware time step. 
For the time of arrival-correlator the computation time is ca. 
14 seconds for an echo integration window of 6 (xbased on a 
1 .7 GHz Intel Xeon processor system. By calculating correla- 
tion coefficients only at the echo peak position the processing 
time decreases to 10 seconds. The FFT correlation function 
contains all available lag times and computation time is in- 
versely proportional to the sampling time t s . For the same 
data it takes 1.7 seconds with 12 /is sampling time and less 
than one second for 24 p.s. Processing the data with the lin- 
ear correlator can be even faster. The whole analysis takes 
approximately 1 .5 seconds if a single channel for each calcu- 
lated echo with a sampling time of 12 [is is choosen. Most 
of this time is needed to convert the raw data to the format 
photon-counts per sampling time, which could be easily done 
on-the-fly during data recording. In conclusion we find that 
the time needed for an optimized data processing scheme is 
negligible compared to the total measurement time. 

Finally we would like to comment on the accessible time 
range at a given rotation frequency. The echo period can not 
be known with absolute accuracy and at very high orders the 
sampling time and echo period will not match any more and 
the signal is lost. We commonly perform measurements up to 
echo number 1000 covering three orders of magnitude in lag- 
time. Since the echo width is only of the order of 1 (xs the pe- 
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FIG. 3: Normalized ICF in backscattering (VH- geometry: perpen- 
dicular polarization). Solid line - ICF from time averaging over 20 
minutes, Symbols - echo analysis of a 12 second measurement (o) - 
data from time-of-arrival data processing, (a) - linear correlator with 
sampling time 12 [is. Inset: Lin-log plot of the same data 

riod has to be known with nanosecond accuracy. To overcome 
this difficulty for even higher order echos one might think of 
monitoring the echo period continuously during processing. 
In this case, however, the echo shape has to be resolved in de- 



tail which goes at the expense of computation speed. A more 
practical way to increase the time window is by performing 
several measurements at different rotation speed. For such a 
scheme the integral time of measurement will still be set ap- 
proximately by the longest measurement. 

In the summary we have shown that our Two-Cell Echo ap- 
proach allows to measure the ensemble averaged DWS corre- 
lation function nearly in real time. Existing DWS experiments 
can be easily upgraded if the laser power is sufficient to drive 
the two-cell echo experiment. Besides a simple device for pre- 
cise mechanical oscillation or rotation and a suitable photon 
counter or correlator no other hardware is needed. Further- 
more, combined with traditional two-cell DWS at very low 
rotation speeds, the full intensity autocorrelation function can 
be measured with a single experimental setup covering more 
than 10 decades in correlation time. 
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